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THE SHAPES OF MILKY WAY SATELLITES: LOOKING FOR SIGNATURES OF TIDAL STIRRING 
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ABSTRACT 

We study the shapes of Milky Way satellites in the context of the tidal stirring scenario for the 
formation of dwarf spheroidal galaxies. The standard procedures used to measure shapes involve 
smoothing and binning of data and thus may not be sufficient to detect structural properties like 
bars, which are usually subtle in low surface brightness systems. Taking advantage of the fact that 
in nearby dwarfs photometry of individual stars is available we introduce discrete measures of shape 
based on the two-dimensional inertia tensor and the Fourier bar mode. We apply these measures of 
shape first to a variety of simulated dwarf galaxies formed via tidal stirring of disks embedded in 
dark matter halos and orbiting the Milky Way. In addition to strong mass loss and randomization of 
stellar orbits, the disks undergo morphological transformation that typically involves the formation 
of a triaxial bar after the first pericenter passage. These tidally induced bars persist for a few Gyr 
before being shortened towards a more spherical shape if the tidal force is strong enough. We test this 
prediction by measuring in a similar way the shape of nearby dwarf galaxies, satellites of the Milky 
Way. We detect inner bars in Ursa Minor, Sagittarius, LMC and possibly Carina. In addition, six 
out of eleven studied dwarfs show elongated stellar distributions in the outer parts that may signify 
transition to tidal tails. We thus find the shapes of Milky Way satellites to be consistent with the 
predictions of the tidal stirring model. 

Subject headings: galaxies: dwarf - galaxies: Local Group - galaxies: fundamental parameters - 
galaxies: kinematics and dynamics - galaxies: structure 



1. INTRODUCTION 

The dwarf galaxies of the Local Group (for reviews see 
Mateo 1998, Tolstoy et al. 2009) display a variety of 
shapes. While this is particularly true for dwarf irreg- 
ular (dlrr) galaxies, it also applies to dwarf spheroidals 
(dSph). Unlike globular clusters, they tend to exhibit 
significantly non-spherical images of stellar content. The 
origin of these shapes is a subject of ongoing debate but 
one of the interesting possibilities involves the evolution- 
ary connection between dlrrs and dSphs. This scenario 
proposes that dSphs evolved from the primordial dlrrs 
via tidal stirring caused by the presence of a large host 
galaxy, such as the Milky Way or Andromeda (Mayer et 
al. 2001). 

The general picture painted by the tidal stirring model 
studied by the means of high resolution W-body simu- 
lations (Mayer et al. 2001; Klimentowski et al. 2009; 
Kazantzidis et al. 2011) involves strong mass loss and 
morphological transformation, as well as a rebuilding of 
the orbital structure of the initial disk. In such simula- 
tions, a primordial disk galaxy embedded in an extended 
dark matter halo is placed on different orbits around a 

1 Nicolaus Copernicus Astronomical Center, 00-716 Warsaw, 
Poland; lokas@camk.edu.pl 

2 Department of Astronomy, University of Virginia, Char- 
lottesville, VA 22904-4325, USA 

3 Center for Cosmology and Astro-Particle Physics; and De- 
partment of Physics; and Department of Astronomy, The Ohio 
State University, Columbus, OH 43210, USA 

4 Institute for Theoretical Physics, University of Zurich, CH- 
8057 Zurich, Switzerland 

5 Department of Physics, Applied Physics, & Astronomy, 
Rensselaer Polytechnic Institute, Troy, NY 12180-3590, USA 

6 Jet Propulsion Laboratory, California Institute of Technol- 
ogy, Pasadena, CA 91109, USA 



Milky Way-sized galaxy and evolved for a significant frac- 
tion of the Hubble time. A thorough study of different 
aspects of the tidal stirring scenario, for different ini- 
tial dwarf structures and different orbits, was recently 
performed by Kazantzidis et al. (2011) and Lokas et al. 
(2011). An important and general feature of the morpho- 
logical evolution of the dwarf established in these stud- 
ies is the formation and longevity of a triaxial, tidally 
induced bar. The bar, usually formed out of the disk 
at first pericenter passage, tumbles with a period much 
shorter than the orbital period of the dwarf, which leads 
to a number of observational biases in the characteris- 
tic parameters of the dwarfs (Lokas et al. 2011). As 
the tidal evolution progresses, the bar shortens and the 
stellar component of the dwarf becomes more and more 
spherical. If the tidal force is strong enough (on tight 
orbits or orbits with small pericenters) , the final shape 
of the stellar component can be almost exactly spherical. 

In this paper we discuss selected observational char- 
acteristics of the tidal stirring model that may help to 
distinguish it from other scenarios for the formation of 
dSph galaxies, e.g., those in which they form in isolation 
by purely baryonic processes such as cooling, star forma- 
tion, feedback from supernovae and UV background radi- 
ation (Ricotti & Gnedin 2005; Tassis et al. 2008; Sawala 
et al. 2010). If the overall picture suggested by the tidal 
stirring model is correct, some of the less evolved dwarfs 
in the Local Group may still be in the bar-like phase. 
Such elongated shapes are indeed observed at least in a 
few cases, such as Ursa Minor (e.g., Irwin & Hatzidim- 
itriou 1995), Sagittarius (Majewski et al. 2003) and the 
more recently discovered Hercules dwarf (Coleman et al. 
2007). 
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TABLE l 

Properties of the simulated dwarfs. 
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Fig. 1. — The effect of smoothing on the detection of the bar-like shapes in dSph galaxies. The left panel shows 2000 stars randomly 
selected from all stars within r < 2 kpc from the center of the dwarf at the apocenter that occurred at 4.2 Gyr from the start of simulation 
Sll when the shape is most bar-like. The positions of the stars are projected on the surface of the sky assuming that the dwarf is viewed 
along the intermediate (y) axis, i.e., perpendicular to the bar. The remaining panels (from the left to the right) show the surface density 
of the stars measured in bins of 0.4 kpc on a side after application of Gaussian smoothing with standard deviation of a = 0.1,0.2 and 0.3 
kpc respectively. The contours are plotted at constant numbers of stars TV with a step of log TV = 0.2 starting with log TV = 2, 1.8 and 1.6 
at the innermost contour for a = 0.1, 0.2 and 0.3 respectively. 



The number of such detections may be low however 
for at least four reasons. First, the bars are oriented 
randomly with respect to an observer located within the 
Milky Way, so some bar-like objects may appear as only 
slightly flattened (which is indeed the case, given that 
the average ellipticity of dSph galaxies is about 0.3). Sec- 
ond, the bars may escape detection due to the smooth- 
ing procedure usually applied when measuring the sur- 
face density distribution of stars in these objects. Third, 
most dwarf galaxies are faint (thus containing a low num- 
ber of stars to sample) and observed against the fore- 
ground of Milky Way stars, which can be only partially 
removed using typical cleaning procedures, such as color- 
magnitude or color-color cuts. Fourth, the intrinsic bar- 
like shapes may be difficult to distinguish from the tidal 
extensions in the surface distribution of the stars ex- 
pected at close encounters between the dwarfs and the 
larger host galaxy. 

The measurements of the shapes of dwarf galaxies in 



the vicinity of the Milky Way are usually performed us- 
ing standard tools that involve binning and smoothing of 
the stellar distribution. In this work we try to improve 
these procedures taking advantage of the fact that the 
satellites of the Milky Way are actually the best extra- 
galactic systems where individual stars can be detected 
and their spatial distribution studied. We introduce dis- 
crete measures of the shape involving only the positions 
of the stars on the surface of the sky. We first apply 
these procedures to simulated dwarfs to see what shapes 
are expected within the tidal stirring scenario. When ap- 
plying the procedures to real satellites of the Milky Way 
we use the cleanest possible stellar samples, selected us- 
ing the information provided by color-magnitude as well 
as color-color diagrams for the stellar populations. 

The paper is organized as follows. In section 2 we 
briefly describe the simulations used in this study. In 
section 3 we present our motivation for using the dis- 
crete measures of shape and introduce them. Section 4 
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Fig. 2. — The shapes of the stellar component in simulation Ol as 
a function of time. The upper panel shows the axis ratios b/a (solid 
line), c/a (dashed line) and c/b (dotted line), which correspond to 
the line-of-sight axis ratio B/A when the dwarf is observed along 
the 2, y and x axis respectively. The lower panel shows the mea- 
surements of the bar parameter A2 when the dwarf is seen along 
the 2 axis (solid line) y axis (dashed line) and x axis (dotted line) as 
a function of time for simulation Ol. Vertical dotted lines indicate 
pericenter passages. The principal axes and the shape values were 
determined for stars within 2r 1 / 2 , where rj/2 is the 3D half-light 
radius determined at the end of the simulation. The final values of 
the shape parameters are listed in Table [T] 

presents the application of these measures to the simu- 
lated dwarf galaxies formed via tidal stirring for a vari- 
ety of initial configurations. In section 5 we measure the 
shapes of the real dwarf galaxies found in the vicinity 
of the Milky Way in an analogous way and compare the 
results to predictions from simulations. The discussion 
follows in section 6. 

2. THE SIMULATIONS 

In this work we use the set of collisionless TV-body sim- 
ulations presented in detail by Kazantzidis et al. (2011) 
and Lokas et al. (2011). The purpose of these stud- 
ies was to elucidate the formation of dSph galaxies via 
tidal interactions between late-type, rotationally sup- 
ported dwarfs and Milky Way-sized hosts. The simu- 
lations, together with all the parameters relevant for the 
present study, are listed in Table [T] with the notation fol- 
lowing that of Lokas et al. (2011). The simulations were 
performed using the multistepping, parallel, tree iV-body 
code PKDGRAV (Stadel 2001). 

Numerical realizations of dwarf galaxy models were 
constructed using the method of Widrow & Dubinski 
(2005). The dwarfs were composed of exponential stel- 
lar disks embedded in cuspy, cosmologically motivated 
Navarro et al. (1996, hereafter NFW) dark matter ha- 
los. The default dwarf galaxy model had a virial mass 
of Mh = 10 9 M Q and a concentration parameter c = 20. 
The disk mass fraction, md, and the halo spin parame- 
ter, A, were equal to 0.02 and 0.04, respectively. The 
resulting disk radial scale length was Rd = 0.41 kpc 
(Mo et al. 1998) and the disk thickness was specified by 
the thickness parameter z^/Rd = 0.2, where Zd denotes 
the (sech 2 ) vertical scale height of the disk. The ori- 



entation of the internal angular momentum of the dwarf 
with respect to the orbital angular momentum was mildly 
prograde and equal to i — 45°. Each dwarf galaxy 
model contained = 10 6 dark matter particles and 
Nd = 1.2 x 10 6 disk particles. The gravitational soften- 
ing was set to e n = 60 pc and £d = 15 pc for the particles 
in the two components, respectively. 

The Milky Way model was constructed as a live real- 
ization of the MWb model in Widrow & Dubinski (2005), 
which satisfies a broad range of observational constraints 
for the present day Galaxy and consists of an exponential 
stellar disk of mass M D = 3.53 x 10 10 M Q , a Hernquist 
(1990) bulge of mass M B = 1.18 xlO 10 M©, and an NFW 
dark matter halo of mass M H = 7.35 x 10 11 M . The 
TV-body realization of the Milky Way used = 10 6 
particles in the disk, Nb = 5 x 10 5 in the bulge, and 
TVh = 2 x 10 6 in the dark matter halo and employed a 
gravitational softening of £d = 50 pc, £b = 50 pc, and 
en = 2 kpc, respectively. 

The default dwarf galaxy model was placed on seven 
different orbits (01-07 in Table [lj around the primary 
galaxy with orbital apocenters, r apo , and pericentcrs, 
r por i, listed in the third and fourth column of Table [TJ In 
simulations S6-S17 (see Table [1]), we varied the orienta- 
tion of the dwarf galaxy disk and the structural parame- 
ters of the dwarfs while keeping them on the same orbit 
(that used for run Ol, with r apo = 125 kpc and r pcrl = 25 
kpc). The second column of Table [T] lists the parameters 
that were varied in each case (in parentheses we give the 
value by which a given parameter was changed). In all 
simulations, the dwarfs were initially placed at the apoc- 
enters of their orbits and the evolution was followed for 
10 Gyr. The orbital times corresponding to these orbits 
are listed in the fifth column of the Table. 

3. MEASURING SHAPES 

The usual approaches to measuring shapes of the Local 
Group dwarf galaxies involve procedures such as binning 
and smoothing of the stellar positions and fitting ellipses 
to the 2D distribution of stars. The results of such pro- 
cedures are typically quantified in terms of the ellipticity 
(or axis ratio) and major axis orientation as a function of 
radius. For most dwarfs, a single ellipticity parameter, 
e = l~b/a (with a and b the major and minor axis length 
respectively) is measured (e.g., Mateo 1998), primarily in 
the outer parts of the dwarf. When the number of stars 
available is low, their positions are often smoothed by 
applying a Gaussian smoothing. If the smoothing scale 
is not carefully chosen, such procedures may cause the 
dwarfs to appear more circular than they really are, re- 
moving substructures such as bars, or at least making 
them less pronounced. 

We illustrate this effect in Figure Q] by performing a 
similar procedure on one of our simulated dwarfs. For 
this purpose we selected stars within radius r < 2 kpc 
from the center of the dwarf (so that there is no contami- 
nation from the tidal tails) from the output of simulation 
Sll when the dwarf is at apocenter (after 4.2 Gyr from 
the start of the simulation) and its shape is most bar-like. 
From this sample we randomly chose 2000 stars whose 
positions (plotted as dots in the left panel of Figure [1} 
were projected on a plane perpendicular to the line of 
sight. The line of sight was chosen to be perpendicu- 
lar to the longest and shortest axis of the bar where the 



4 



E. L. Lokas et al. 



0.5 











1. 

0.5 






02 














0. 



0.5 



0.5 



0.5 



0.5 



0.5 



















( 


35 



























S6 



















S9 



















S12 























si; 





2 4 6 

time [Gyr] 







'03 ; 







10 



10 



8 10 



0.5 



10 



0.5 - 



0. 

10 



















04 



J 0. 
10 

























sio ; 



0.5 



10 




2 4 6 

time [Gyr] 



10 











07 




































S8 





10 



















Sll 







10 





















S14 





10 

























S17 



2 4 6 

time [Gyr] 



10 



Fig. 3. — The axis ratios b/a (solid line), c/a (dashed line) and c/b (dotted line) for simulations 02-S17. The values correspond to the 
line-of-sight axis ratio B/A when the dwarf is observed along the z, y and x axis respectively. The values were determined for stars within 
2^1/2 with )-j/2 adopted to be the final value for each simulation. Vertical dotted lines indicate pericenter passages. The final values of the 
shape parameters are listed in Table Q] 



detection of the bar should be easiest. 

We then replaced each star by a 2D Gaussian centered 
on the star and normalized to unity with a standard de- 
viation a. The observed region was divided into bins of 
0.4 x 0.4 kpc and the number of stars in each bin was cal- 
culated by summing up the contributions from all Gaus- 
sians. The resulting surface density distributions of the 
stars obtained with a = 0.1,0.2 and 0.3 kpc are shown 
respectively in the second, third and fourth panel of Fig- 
ure Q] We can see that while with a = 0.1 kpc the bar 
is quite well visible, with er = 0.3 kpc the shape would 
hardly be classified as bar-like and more probably de- 
scribed as highly elliptical. Note that the values of the 
smoothing parameter a adopted here are quite realistic, 



i.e., similar to what is usually applied to real data. For 
example, in studies of the Leo I dSph galaxy Sohn et 
al. (2007) and Mateo et al. (2008) used a = 1.1 and 
2 arcmin respectively, which corresponds to 0.2 and 0.4 
kpc, whereas the half-light radius of Leo I is only 0.25 
kpc. We learn from this exercise that when the smooth- 
ing scale becomes comparable to the half-light radius of 
the dwarf (r 1 / 2 = 0.48 kpc for our simulated dwarf at the 
selected time) , interesting morphological features such as 
bars may no longer be detectable even if the line of sight 
is perpendicular to the bar. 

To avoid such problems, in this study we will quan- 
tify shapes using only discrete measures, i.e., ones that 
use the 2D positions of the stars, without applying any 
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Fig. 4. — The measurements of the bar parameter A2 when the dwarf is seen along the z axis (solid line) y axis (dashed line) and x axis 
(dotted line) as a function of time for simulations 02-S17. The values were determined for stars within lT\ji with ri/2 adopted to be the 
final value for each simulation. Vertical dotted lines indicate pericenter passages. The final values of the shape parameters are listed in 
Tabled] 



smoothing or binning procedures. We also aim to mea- 
sure the intrinsic shapes, i.e., in the absence of any con- 
tamination from the Milky Way or tidally stripped stars. 
Although avoiding such contamination is not possible in 
real observations, when talking about our comparison 
shapes from the simulations, we will always measure such 
intrinsic shapes using stars found within 2r 1 / 2 from the 
center of the simulated dwarf, where r\/2 is the 3D half- 
light radius containing half the total number of stars in 
3D at the end of the simulation, as discussed in Lokas et 
al. (2011). The values of these radii for our models are 
listed in the sixth column of Table [TJ These fixed radii 
are used for all outputs and all lines of sight whenever the 
shapes of the simulated dwarfs are discussed. Note that 
such ri/2 values arc different from their 2D counterparts 



that in general depend on the line of sight. 

We will use two measures of the shape, both calculated 
from the 2D distribution of stars within a fixed radius. 
In each case we first determine the principal axes of the 
stellar distribution from the 2D inertia tensor and ro- 
tate the stars so that their x(y) positions lie along the 
major (minor) axis. The first measure of the shape we 
apply is the axis ratio B/A, or alternatively the ellip- 
ticity e = 1 - B/A, where B/A = [Eij/f/EiX?] 1 / 2 . The 
second measure of the shape we use is the so-called bar 
amplitude of the Fourier decomposition of stellar phases, 
A2 = \T,j exp(2icf>j)\/N where <f>j are the phases of the 
positions of stellar particles in the xy plane. As we will 
demonstrate below, these measures, although slightly dif- 
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Fig. 5. — The evolution of the axis ratio c/b versus b/a for dwarfs 
on different orbits 01-07 (upper panel) and dwarfs with different 
initial structure S6-S17 (lower panel). Color lines join measure- 
ments at subsequent apocenters starting from the initial values 
at the lower right region of each panel. The diagonal black line 
marks the transition between the oblate (b/a > c/b) and prolate 
(b/a < c/b) regimes. 

ferent in their numerical values, are highly correlated and 
can be used alternatively. 

4. THE SHAPES OF THE SIMULATED DWARFS 

To measure the shapes of the stellar components of the 
simulated dwarfs, for each output of each run listed in 
the first column of Table [1] we selected stars contained 
within 2r\/2 from the dwarf's center. The 3D half-light 
radius r\/2 was kept at a fixed value for all outputs of 
a given run, and equal to its final value at the end of 
the simulation, i.e., after 10 Gyr of evolution. Although 
r 1 / 2 depends slightly on time (see Lokas et al. 2011), by 
keeping the radius fixed we make sure that the shapes are 
measured at the same scale for all outputs and therefore 
that their evolution in time is easier to interpret (note 
that in Kazantzidis et al. 2011 the shapes for the same 
runs were determined within a variable radius r max cor- 
responding to the maximum circular velocity V max ). 

For each sample of stars selected in this way, we deter- 



mine the orientation of the principal axes from the 3D 
inertia tensor and rotate the distribution to align stellar 
positions with those axes, so that x is along the longest 
axis a, y along the intermediate axis b and z along the 
shortest axis c. The 3D shapes of the dwarfs can then 
be specified by the axis ratios b/a = [Ei^/Eja; 2 ] 1 / 2 , 
c/a = [Eizf/Eixf] 1 / 2 and c/b = [EiZ 2 /E i2 / 2 ]V2 where 
(xi, yi, zi) is the position of the zth star. The 2D shapes 
are then obtained by selecting a line of sight. If the lines 
of sight are along the z, y and x axis, then the measured 
2D axis ratio B/A will be equal to b/a, c/a and c/b re- 
spectively. An example of the shapes measured in this 
way is shown as a function of time in the upper panel 
of Figure for simulation 01. Using the same lines of 
sight we measure the bar mode parameter Ai which is 
shown in the lower panel of this figure. Analogous results 
for the remaining runs 02-07 and S6-S17 are plotted in 
Figures [3] and HI 

The results presented in Figures [2]|4] confirm those ob- 
tained in a slightly different way by Kazantzidis et al. 
(2011): in all cases after the first pericenter the stellar 
component transforms from a disk to a triaxial ellipsoid 
with a =/= b =/= c. In most runs, at some point in the evo- 
lution, the shape becomes prolate with c/b > b/a which 
can be interpreted as the formation of a bar-like shape. 
This is confirmed by high values of A^ > 0.2 detected 
for the line of sight along the z axis, i.e. signifying the 
formation of a bar in the plane of the initial disk. As 
the evolution progresses, the shape becomes more spher- 
ical with all axis ratios approaching unity, with the most 
significant changes occurring at the pericenter passages 
(vertical dotted lines in Figures [2E| when the dwarfs are 
most strongly tidally shocked. 

The evolution of the global shape of the stellar compo- 
nent can be illustrated further by plotting the evolution- 
ary tracks of the dwarfs in the b/a — c/b plane. We show 
them in Figure [5] with lines of different color for each 
simulation (see the last column of Tabled]) for runs 01- 
07 in the upper panel and runs S6-S17 in the lower one. 
The tracks start at the lower right corner of each panel 
and progress towards the upper left in subsequent apoc- 
enters. Most of the tracks cross the diagonal black line at 
some point, signifying the transition from the oblate to 
the prolate shape, which can be interpreted as the forma- 
tion of a bar. These correspond to the dotted line being 
above the solid line in the upper panel of Figure [2] and 
in Figure [3] Interestingly, all transitions into the prolate 
half of the diagram are accompanied by significantly high 
values of A^ measured along the shortest (z) axis (solid 
line) in the lower panel of Figure [2] and in Figure QJ This 
is due to the fact that if the bar mode occurs in the ini- 
tial disk plane usually it means that the global shape of 
the dwarf becomes bar-like, i.e., most stars participate 
in the bar instability and form a bar while very few stars 
remain elsewhere. On the other hand, in cases when the 
shape remains oblate during the whole evolution (runs 
03, 05, 06, 07, S10, S15 and S16) the z-axis bar mode 
remains low with Ai < 0.2. An interesting borderline 
case is S10 when a bar mode is weakly excited (A2 ~ 0.2 
between the second and third pericenter) but the over- 
all shape of the dwarf remains oblate (c/b < b/a at all 
apocenters). 

In addition to the z-axis bar mode shown in the lower 
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Fig. 6. — The surface density distribution of the stars seen in simulation Ol at the second, third, fourth and fifth apocenter (from the 
left to the right column) along the principal axes of the stellar distribution x, y and z (from the upper to the lower row). The bar formed 
after the first pericenter, positioned along the x axis, is gradually diminished at later stages. 

01 the dwarf is seen along the y and x axis of the stellar 

distribution. Obviously, a distant observer will see the 
dwarf along a random line of sight and will not be able to 
determine a full 3D shape of the stellar component. We 
thus show what such measurements would give in terms 
of A2 if the line of sight was close to the y or x axis. 
In these cases, the dwarf may often appear elongated 
not because it is a bar but because it is still a disk, as 
demonstrated by high values of A2 at the early stages of 
evolution. Our observer may not be able to distinguish 
a disk from a bar even if kinematical measurements were 
available because, according to our simulations, a tidally 
induced bar is expected to tumble (see Figure 4 in Lokas 
et al. 2011) and thus may display a similar rotation curve 
as a disk. Only additional precise photometric measure- 
ments aiming to constrain the depth of the dwarf galaxy 
could differentiate between a disk and bar in this case. 

Let us now focus on the 2D shapes of the simulated 
dwarfs, as they would be seen by a distant observer. We 
illustrate the evolution of such shapes using our default 
simulation 01. In Figure [5] we show how the dwarf's 
stellar component would appear when viewed along the 
principal axes of the stellar distribution x, y and z (from 
the upper to the lower row) at subsequent apocenters of 
its orbit around the Milky Way (from the left to the right 
column). The images were produced by selecting stars 
within r < 2 kpc and rotating their positions to align 




x [kpc] 



Fig. 7. — The surface density distribution of the stars seen at the 
final output of simulation Ol along the y (intermediate) axis of the 
stellar distribution, i.e., at its most non-spherical appearance. 

panel of Figure [5] and in Figure 2] with solid lines, we 
also add (with dashed and dotted lines respectively) the 
analogous measurements of this parameter done when 
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Fig. 8. — The surface density distribution of the stars seen at the final output of simulations 02-S17 along the y (intermediate) axis of the 
stellar distribution, i.e., at their most non-spherical appearance. For all simulations the distributions are shown for stars within r < 2 kpc 
although the orientation of the principal axes of the stellar component were determined using stars within r < 2r 1 / 2 with rj/2 different for 
each simulation. The stellar distribution in 04 stands out with its pronounced, elongated outer contours due to the very recent pericenter 
passage. 
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Fig. 9. — The linc-of-sight shape parameters 1 — B/A (filled symbols) and A2 (open symbols) determined from the 2D surface density 
distribution of the stars in the simulated dwarfs. The triangles, squares and circles refer to measurements performed for observations along 
the intermediate (y), longest (x) and shortest (z) axis of the stellar distribution. The plotted values are listed in Table[TJ 



TABLE 2 

Properties of the Milky Way satellites included in this analysis. 



Dwarf 
galaxy 


Center 
a(J2000) 


Center 
5(J2000) 


d Q 
[kpc] 


n/2 

[arcmin] 


n/2 

[kpc] 


-^"stars 

(< 2r 1/2 ) 


Selection 
method 


Source 
of data 


Fornax 


02 ii 39 m 4gs 


-34°31'01" 


138 


16.3 


0.656 


6283 


CM 


1,2,1 


Sculptor 


01 h 00 m 06 s 


-33°42'40" 


86 


8.85 


0.221 


1751 


CM+CC 


3,4,3 


Carina 


06 h 41 m 37 s 


-50°58'00" 


105 


7.86 


0.240 


269 


CM+CC 


5,6,7 


Ursa Minor 


15 h 09 m 04 s 


+67°13'51" 


76 


14.4 


0.318 


465 


CM+CC 


8,9,10 


Sextans 


10 h 13 m 03 s 


-01°36'54" 


95 


27.9 


0.770 


2503 


CM 


5,11,12 


Draco 


17 h 20 m 19 s 


+57°54'48" 


80 


9.66 


0.225 


1466 


CM 


5,13,12 


Leo I 


10 h 08 m 29 s 


+12°18'20" 


257 


3.34 


0.250 


1100 


CM+CC 


14,14,14 


Leo II 


ll h 13 m 29 s 


+22°09'06" 


233 


2.28 


0.155 


272 


CM+CC 


15,16,17 


Sagittarius 


18 h 54 m 59 s 


-30°27'38" 


25 


221 


1.604 


3344 


CM+CC 


18,19,18 


LMC 


05 h 23 m 34 s 


-69°45'24" 


50 


136 


1.979 


132549 


CM 


5,20,21 


SMC 


00 h 51 m 00 s 


-73°07'12" 


62 


70.9 


1.278 


13626 


CM 


22,23,21 



References. — (1) Coleman et al. (2005); (2) Saviane et al. (2000); (3) Westfall et al. (2006); (4) Pietrzyriski et al. (2008); (5) Mateo 
(1998); (6) Pietrzyriski et al. (2009); (7) Majewski et al. (2000b); (8) Kleyna et al. (1998); (9) Bellazzini et al. (2002); (10) Palma ct al. 
(2003); (11) Lee et al. (2003); (12) SDSS DR8, Aihara et al. (2011); (13) Aparicio et al. (2001); (14) Sohn et al. (2007); (15) Coleman et 
al. (2007); (16) Bellazzini et al. (2005); (17) Siegel & Majewski (2000); (18) Majewski et al. (2003); (19) Kunder & Chaboyer (2009); (20) 
Schaefer (2008); (21) 2MASS survey, Skrutskie et al. (2006); (22) Gonidakis et al. (2009); (23) Szewczyk et al. (2009). 



TABLE 3 

The shape parameters of Milky' Way satellites included in this analysis. 



Dwarf 




1 - B/A 






A 2 




galaxy 


whole 


inner 


outer 


whole 


inner 


outer 


Fornax 


0.18 


0.09 


0.21 


0.12 


0.07 


0.22 


Sculptor 


0.14 


0.07 


0.17 


0.11 


0.06 


0.18 


Carina 


0.20 


0.17 


0.21 


0.20 


0.16 


0.24 


Ursa Minor 


0.38 


0.34 


0.39 


0.36 


0.29 


0.46 


Sextans 


0.07 


0.09 


0.07 


0.08 


0.08 


0.08 


Draco 


0.16 


0.15 


0.16 


0.15 


0.14 


0.18 


Leo I 


0.23 


0.08 


0.27 


0.14 


0.04 


0.26 


Leo II 


0.09 


0.10 


0.09 


0.09 


0.10 


0.09 


Sagittarius 


0.47 


0.32 


0.51 


0.39 


0.26 


0.61 


LMC 


0.07 


0.19 


0.04 


0.15 


0.20 


0.08 


SMC 


0.18 


0.13 


0.19 


0.14 


0.10 


0.21 



Note. — For each dwarf galaxy we list three values of parameters 1 — B/A and A 2 : for the whole stellar sample within < r < 2ri/ 2 , 
for the inner sample within < r < t^m and for the outer sample within r-iy 2 < T < 2rj/ 2 . The values of A 2 > 0.2 signifying the presence 
of bar-like shape were boldfaced. 
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them with the principal axes determined for stars with 
r < 2ri/2, as described above. The plots show the surface 
density distribution of the stars with contours equally 
spaced in log density. In agreement with the conclusions 
from Figure [5J at second apocenter the dwarf already 
has a pronounced bar and its shape changes from the 
oblate (initial disk) to the prolate (a bar). At the third 
apocenter the shape is again oblate (as shown by the 
green line in the upper panel of Figure [5] which crosses the 
diagonal line back again from the prolate to the oblate 
part of the diagram at this time). Later on, the shape 
gradually becomes more and more spherical with all axis 
ratios approaching unity. The final image of the dwarf, 
at the end of the simulation, at its most non-spherical 
appearance (i.e., viewed along the y axis) is shown in 
Figure [3 

In Figure [8] we show analogous images for the remain- 
ing simulations 02-07 and S6-S17. In all cases stars 
within r < 2 kpc were extracted from the last output of 
the simulation and rotated to align the distribution with 
the principal axes. Only the views along the y axis are 
shown, i.e. we see the dwarfs along the line of sight in 
which they appear most non-spherical. The dwarfs are 
shown at the end of their evolution when most are far 
from their pericenters, and therefore weakly affected by 
tidal forces at this time, except for the case of 04 which 
has just passed its 6th pericenter. Although this dwarf 
has already reached an almost spherical shape during its 
earlier evolution, it is strongly affected by tidal forces 
from the Milky Way at this particular instant (note that 
its pericentric distance was 12 kpc at this passage) which 
manifests itself by strong distortion of the outer contours. 

Going back to our single-number measures of the 
shapes, 1 — B/A and A2, we summarize our findings by 
comparing in Figure [9] the values of these parameters 
measured along the three lines of sight at the end of each 
simulation and listed in columns 7-12 of Table [TJ For 
each run we plot 1 — B/A with filled symbols and Ai with 
open symbols. The triangles, squares and circles corre- 
spond to measurements along the y, x and z axis of the 
stellar component, respectively. As discussed above, the 
observation along the intermediate axis y always gives 
the most non-spherical appearance, i.e., triangles always 
mark the highest values of 1 — B/A and A2. In addition, 
the next highest value of cllipticity is almost always the 
one seen along the x axis (so that the points are ordered 
triangle-square-circle from top to bottom). The excep- 
tions occur for runs 04, S7, Sll and S12. Since these 
cases correspond to b/a < c/b this means that only in 
these four cases do the dwarfs end up more prolate than 
oblate (04 should be excluded because its shape is due 
to strong tidal forces at pericenter). It is interesting, 
therefore, that although almost all dwarfs go through a 
bar-like phase during their evolution, most of them end 
up more oblate than prolate. 

We also note that the values of 1 — B/A and A2 are 
quite similar, although not exactly equal, and they pre- 
serve the relative triangle-square-circle hierarchy of the 
measurements. To illustrate this similarity further, in 
Figure fTOl we plot 1 — B/A as a function of A2 for the 
line of sight where the dwarfs appear most non-spherical 
(observation along y) . The lines join the values measured 
at subsequent apocenters, from the most non-spherical at 
the upper right corner of each panel to the almost spher- 
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Fig. 10. — The correlation between different measures of shape, 
1 — B/A and A2, for the line of sight along the y axis (when the 
dwarfs appear most non-spherical) for dwarfs on different orbits 
01-07 (upper panel) and dwarfs with different initial structure 
S6-S17 (lower panel). Color lines join measurements at subsequent 
apocenters starting from the initial values at the upper right region 
of each panel. 

ical in the lower left. The upper panel shows the results 
for simulations of dwarfs on different orbits (01-07) and 
the lower one those with different structure of the dwarf 
(S6-S17). The lines corresponding to each simulation can 
be identified by referring to the last column of Table [TJ 
We can see that there is a one-to-one correspondence 
between the two measures of shape, and thus conclude 
that both measures, although defined differently from a 
mathematical point of view, essentially measure the de- 
partures from sphericity in a similar way and can be used 
interchangeably to quantify the shapes of galaxies with 
resolved stellar populations. 

5. THE SHAPES OF MILKY WAY SATELLITES 

In this section we discuss the shape parameters of the 
real dwarfs, satellites of the Milky Way using photo- 
metric data. Our dwarf galaxy sample includes objects 
within heliocentric distance of about 260 kpc (close to the 
largest apocenter considered in the simulations) , brighter 
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than My = — 8 mag and with good photometric data 
available. These requirements restrict the sample to the 
eight classical dSph galaxies, the Sagittarius dwarf and 
the Large and Small Magellanic Clouds (LMC and SMC). 
The data sets we used comprise the positions and mag- 
nitudes of stars in at least two bands to allow for the 
selection of probable members of dwarf galaxies from the 
color- magnitude diagrams (CMD). In the case of the For- 
nax dwarf we used the data from Coleman et al. (2005), 
for Sextans and Draco the data were retrieved from the 
Sloan Digital Sky Survey database. For these dwarfs we 
selected member stars by iterative fitting of the position 
of the red giant branch (RGB) and choosing stars within 
a strip of variable width along the branch whose shape 
was adjusted to match the equal-density contours of the 
stars in the CMD. For the LMC and SMC, we extracted 
the stars from the 2MASS survey and also included stars 
from the asymptotic giant branch (see van de Marel & 
Cioni 2001; Gonidakis et al. 2009). 

For the remaining dwarfs (Sculptor, Carina, Ursa Mi- 
nor, Leo I, Leo II and Sagittarius) we used samples 
obtained previously by S.R.M. and collaborators with 
more precise photometric membership determination us- 
ing color-magnitude as well as color-color diagrams to 
extract only the giant stars. The technique, described in 
detail in Majewski et al. (2000a), relies on the sensitiv- 
ity of the DDOhl filter to magnesium features that re- 
flect stellar surface gravity, temperature and abundance 
in later type stars. When combined with the wide-band 
M and T 2 filters of the Washington system, the DDOhl 
filter is especially useful for discriminating giant stars 
from foreground dwarfs on the basis of differences in their 
respective M — DDObl colors at a given M — T 2 color. 
The data, although not publicly available, are described 
in the references provided in Table [5] and we refer the 
reader there for details concerning each dwarf galaxy. 

Table [5] lists the properties of the dwarfs included in 
this analysis: the dwarf name (column 1), the position 
of the adopted center (columns 2 and 3), the heliocentric 
distance of the dwarf (column 4), the 2D half-light radius 
r-yji in arcmin and kpc (columns 5 and 6), the number 
of member stars within 2r 1 / 2 (column 7), the method 
used for the membership determination (column 8, CM = 
color-magnitude, CM+CC = color-magnitude plus color- 
color), and the source of data (column 9, the three num- 
bers give respectively the reference for the dwarf center, 
distance and the photometric data used for the measure- 
ments of its shape). 

Although the values of the half-light radii are avail- 
able in the literature, for consistency we determined them 
from the surface density distribution of the selected mem- 
ber stars by fitting a Plummer profile. For this purpose, 
the surface density profiles were measured in logarith- 
mic bins of projected radii and fitted with data points 
weighted by density adjusting two parameters, the half- 
light radius and normalization. The obtained values of 
ri/2, listed in Tabled are in good agreement with pub- 
lished scales (see, e.g., Table 2 in Lokas et al. 2011). 

We measured the shape parameters 1 — B/A and A 2 
first using member stars in the whole range < r < 
2r 1 / 2 - The values of these parameters are listed in the 
second and fifth column of Table [3] respectively and 
marked as 'whole'. In order to distinguish between intrin- 



sic bars and external tidal extensions it is also interesting 
to look at the parameters calculated separately in the in- 
ner part within < r < 7*1/2 and the outer part within 
r i/2 < r < 2r 1 / 2 - Those are shown in the columns of Ta- 
ble [3] marked respectively as 'inner' and 'outer'. In the 
following we will focus on A 2 , which is a more customary 
quantity used to describe bars, and adopt A 2 > 0.2 as 
a threshold signifying the presence of a bar-like shape. 
Values above this threshold were typed in boldface in 
Tabled 

The motivation behind this division into inner and 
outer sample is that the tidal effects on the shape of 
dwarf galaxies can take two forms. The shapes can be 
affected by tides so that the whole orbital structure of 
the dwarf is rebuilt and a bar forms, or only the outer 
contours of the stellar distribution are elongated signify- 
ing the transition to the tidal tails. The latter situation 
is illustrated very well by the special case of simulation 
04 which ends at the pericenter of the orbit. This mo- 
ment of the evolution is shown in the upper right panel 
of Figure [8] While the inner contours of the stellar den- 
sity distribution are circular, the outer ones are strongly 
elongated due to strong tidal force at the pericenter of 
the orbit. On this orbit with a very small pericenter (as 
well as for the very tight orbit 02) the dwarf transforms 
into an almost perfectly spherical shape early on in the 
evolution, as demonstrated by the time-dependence of 
the shape parameters in Figures [3JU The shape remains 
spherical for most of the time at later stages, except at 
pericenters, where strong elongation takes place. 

In Figure [TTJ we show the contour maps of the surface 
density of the stars in real dwarfs up to projected radii of 
the order of 2 7-1/2- The maps were obtained by smoothing 
the discrete distribution of stars with Gaussians of scale 
a = 0.2ri/2 and binning it into 10 x 10 grids. As discussed 
in section 3, such smoothing scales are small enough not 
to obscure subtle morphological features. The contours 
joining the regions of equal surface density of the stars 
are equally spaced in log surface density. 

As verified in Table 02 only in the case of Ursa Mi- 
nor and Sagittarius are all three values of the bar mode 
parameter — for the whole, inner and outer samples — 
found to have A 2 > 0.2. This strong elongation, deter- 
mined by the discrete shape measures, is confirmed in 
a more qualitative way by the images of these dwarfs 
in Figure [TT] We conclude that these dwarfs possess the 
tidally induced bar and that their outer parts are affected 
by tides. In the case of the Sagittarius dwarf, the outer 
parts are likely stretched by tidal forces acting now, since 
the dwarf is very close to the pericenter of the orbit and 
the pericenter distance is very small (< 15 kpc); (like 
the dwarf in simulation 04 at the end of the evolution). 
In the case of Ursa Minor, which is much farther away 
from the center of the Milky Way at present, the outer 
contours probably signify a transition to the tidal tails. 
For the remaining dwarfs, the inner bar is detected only 
in the LMC and possibly in Carina (A 2 — 0.16 in the 
inner part). In addition, Fornax, Carina, Leo I and the 
SMC have A 2 > 0.2 in the outer parts, signifying mild 
tidal extensions. 

6. DISCUSSION 

We have studied the shapes of dwarf galaxies orbiting 
the Milky Way using new discrete measures appropriate 
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Fig. 11. — The surface density distribution of stars in Milky Way satellites. In all cases (except for the LMC and SMC, where the number 
of stars is large) the surface density of stars was smoothed with Gaussian filters of scale 0.2r]y 2 (where T\ji is the half-light radius from 
Table [2j and the maps were created by binning the resulting stellar distribution into a 10 X 10 mesh. The contours are equally spaced in 
log surface density distribution. The images were rotated to align the major axis (as determined within 2r 1 / 2 ) with the x axis of the plots. 
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for systems where individual stars are resolved. First, 
we traced the evolution of shape parameters in a set of 
collisionless TV-body simulations of late type progenitors 
of present day dSph galaxies. The initially disky dwarfs 
were placed on different orbits around a Milky Way-type 
galaxy and the inner structure of the satellites was also 
varied. In addition to the strong mass loss and ran- 
domization of stellar orbits, tidal interactions inherent 
to such configurations lead to significant evolution of the 
shape. Usually at the first pericenter passage the dwarfs 
undergo bar instability and their overall distribution of 
stars transforms from the initial disks to triaxial bars. 
The bars are preserved for a few Gyrs of evolution but 
become shorter in time and produce, in the most tidally 
affected cases, almost spherical stellar components. 

The effect of tidal forces on the shape of the dwarfs can 
be twofold. In addition to the tidally induced bar, which 
can be easily identified with the full 3D information avail- 
able in the simulations, the tidal forces can stretch the 
outer contours of the stellar distribution even in dwarfs 
that have already evolved into spheres. This is partic- 
ularly well visible at pericenters but can also occur at 
other parts of the orbit where the dwarf is not being 
strongly shocked, but is already surrounded by the ma- 
terial stripped from its main body in the form of tidal 
tails. The two effects can be distinguished by measur- 
ing shapes separately in the inner and outer parts of the 
dwarfs. 

We found the bar phase of the evolution to be quite 
common among the simulated dwarfs. Given that the 
dwarf models were not susceptible to bar instability in 
isolation, the presence of the bars can be considered as 
a signature of tidal stirring. One has to keep in mind 
however, that the tidal forces are extremely effective not 
only in inducing the bars but also in destroying them. If 
the tidal force is very strong, the bar can transform into 
a spherical shape over a rather short timescale. It would 
therefore be very difficult to find a strong correlation 
between the formation of a bar and the strength of the 
tidal force or a dwarf's distance from the Milky Way. 

Nevertheless, given the variety of shapes within the 
dwarf galaxy population in the vicinity of the Milky Way 
we attempted to check if the tidal stirring scenario is 
broadly consistent with the properties of real dwarfs. To 
this end we measured the shapes of brighter dwarf galax- 
ies having good enough available photometry for large 
samples of stars. Out of eleven dwarfs studied, we found 
clear signatures of inner bars in Ursa Minor, Sagittarius, 
LMC and possibly Carina. In six out of eleven dwarfs we 
found indications of tidal extensions in the outer parts 
which likely signify transitions to the tidal tails. 

The elongated shapes of Ursa Minor, Sagittarius, Ca- 
rina and the inner bar of the LMC have been pointed out 
by many previous studies. The results are however par- 
ticularly interesting for Sagittarius because its elongated 
shape is usually interpreted as due to tidal forces acting 
at the present time when the galaxy is at the pericenter 
of its orbit. We find instead that the shape is elongated 
down to the very center of the dwarf thus lending support 
to the model proposed by Lokas et al. (2010a). In this 
scenario, the progenitor of Sagittarius was a disk galaxy 
similar to the present LMC, which transformed into a bar 
at the first pericenter passage. The bar is preserved until 
the present and the dwarf is currently at the second peri- 



center on its orbit around the Milky Way. The presence 
of the inner bar in the LMC also seems consistent with 
the tidal stirring scenario which predicts that bars are 
typically formed at the first pericenter. Since the LMC 
is probably on its first passage around the Milky Way 
(Bcsla et al. 2007), it should be forming its bar right 
now. 

It may seem surprising that among the eight classical 
dSph galaxies in the vicinity of the Milky Way a clear 
inner bar is detected only in Ursa Minor which, accord- 
ing to the scenario presented here, would imply that it 
has not passed many pericenters on its orbit around the 
Milky Way. On the other hand, Ursa Minor is domi- 
nated by the most metal-poor stellar population among 
the classical Milky Way dSphs (Kirby et al. 2011). This 
means that it lost its gas rather early, probably via ram 
pressure stripping on its first passage around the Milky 
Way (Mayer et al. 2007), a mechanism acting most ef- 
fectively for dwarfs accreted early. However, the gas may 
also have been lost via a strong blow-out in the formation 
of the initial stellar population while the system was still 
evolving in isolation. Given the stochasticity in the star 
formation history outcomes, as well as the wide distribu- 
tion of orbital and structural parameters of the progeni- 
tor dwarfs, it would be difficult to exactly reproduce the 
history of any given dwarf just from its present shape. 
The case of Ursa Minor however proves that bars are 
present even among the eight classical dwarfs and may 
be present also in other systems but not detectable be- 
cause of their orientation along our line of sight. 

The present study was based on photometric data 
alone. However, considerable progress towards discrimi- 
nating between different scenarios of dSphs formation is 
expected if these data are combined with kinematic mea- 
surements. Although rich samples of radial velocity mea- 
surements are presently available only for a few dSphs 
(e.g., Walker et al. 2009), such samples are bound to 
grow substantially in the near future thanks to a new gen- 
eration of multi-object spectrographs (such as the Prime 
Focus Spectrograph on the Subaru telescope). Kinemat- 
ical samples for a few thousand stars could help detect 
remnant rotation that should be present in at least some 
dSphs if they formed from disky progenitors, as envi- 
sioned by the tidal stirring model. Such a rotation signal 
has until now only been tentatively detected for some 
dwarfs. In addition, such data will also help distinguish 
bars from other structures. Bars seen perpendicular to 
the longest axis should show little rotation in the inner 
parts, contrary to disks seen edge on. On the other hand, 
rotation detected in systems that appear spherical may 
indicate a bar seen along the longest axis (see Lokas et 
al. 2010b). 

We end with a note of caution that the picture estab- 
lished here with the help of collisionless simulations using 
standard assumptions concerning, e.g., the initial dark 
matter profiles of the dwarfs, may well be modified if hy- 
drodynamical processes are included or the assumptions 
changed. Although few hydrodynamical simulations of 
tidal stirring have been performed to date, they indicate 
that the shapes may be significantly affected by the pres- 
ence of the gas and star formation (Mayer et al. 2007). 
The details of the morphological transformation will also 
be different if the initial dark matter profile has a core 
rather than a cusp (Governato et al. 2010; Lokas et al. 
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2012) and the initial stellar and gaseous disk of the dwarf 
is thicker. We plan to address these issues in our forth- 
coming papers. 
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